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1.  INTRODUCTION 

In  [5]  R.  C.  Grinold  proposed  a linear  programming  cohort  model  which 
could  be  used  for  studying  the  design  of  a manpower  system.  In  the  present 
paper  it  is  suggested  that  his  linear  programming  model  can  be  approximated 
by  a network  transshipment  model,  in  the  same  way  that  A.  Charnes , W.  W.  Cooper, 

K.  A.  Lewis,  and  R.  J.  Niehaus  suggest  approximating  a certain  goal  programming 
model  with  an  imbedded  Markov  chain  by  a network  model,  see  [1,2].  The  loss 
of  generality  in  such  approximations  is  slight,  while  there  are  substantial 
gains  possible  in  terms  of  computing  speeds  as  shown  by  F.  Glover,  K.  Karney, 
and  D.  Klingman  in  [4]  and  V.  Srinivasan  and  G.  Thompson  in  [10]. 

Besides  the  increase  in  speed  of  computation,  many  other  advantages 
become  evident.  One  is  the  availability  of  the  operator  theory  of  parametric 
programming  for  such  problems  developed  by  Srinivasan  and  Thompson  in  [9] 
which  makes  the  process  of  approximating  a linear  program  by  a network  model 
not  difficult.  It  also  makes  the  computation  of  the  pareto  optimal  effectiveness- 
cost  tradeoff  curve  discussed  in  Section  5 very  easy,  see  [11], 

In  Section  2 the  original  model  by  Grinold  is  briefly  described,  and 
then  translated  into  a network  transshipment  model.  Two  objective  functions, 
minimum  cost  and  maximum  effectiveness  are  considered.  In  Section  3 a simple 
numerical  example  with  minimum  cost  objective  is  solved.  Section  4 discusses 
the  interpretations  of  the  dual  matrix  for  the  example  and  for  the  maximum 
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ef fectiveness  objective.  In  Section  5 the  model  is  considered  with  the 
two  objectives  simultaneously  evaluated,  and  the  problem  of  finding  the 
effectiveness-cost  tradeoff  curve  is  solved.  Section  6 is  devoted  to  career 
analysis  that  is  available  from  the  outputs  to  the  model,  and  Section  7 gives 
conclusions . 

It  is  with  the  greatest  pleasure  that  I dedicate  this  paper  to 
Professor  William  W.  Cooper.  He  was  the  one  who  first  introduced  me  to  the 
manpower  planning  area  and  who  has  consistently  supported  my  work  as  well  as 
that  of  many  others  in  the  area.  His  unselfish  devotion  to  scientific  re- 
search, of  all  kinds,  has  greatly  enriched  every  area  which  he  has  touched. 

2 . THE  MODEL 

We  begin  by  introducing  some  notation  very  similar  to  that  employed 
by  Grinold  in  [5].  After  a brief  discussion  of  the  problem  in  that  notation 
we  will  change  notation  to  make  possible  a network  formulation. 

Consider  the  manpower  system  which  has  ranks  r e {l,2,...,R},  with 
increasing  rank  associated  with  increasing  rank  number,  and  years  of  service 
t e {l,2,...,T}.  Thus  T is  the  maximum  period  of  service  and  retirement 
must  take  place  on  or  before  the  beginning  of  year  T + 1.  We  define  the 
following  quantities: 

1 (r,t)  * state  of  being  in  rank  r serving  in  year  t. 

K - number  of  states  (r,t)  for  1 < r < R and  1 < t < T. 

n(r,t)  ■ number  of  people  in  rank  r serving  in  year  t. 

There  are  two  terminal  states,  Retired  and  Separated. 
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f [ (r,t); (r ,t+l)]  * fraction  of  people  in  state  (r,t)  going  to  state  (r,t+l) 
f [ (r , t);  (r+1 , t+1)]  = fraction  of  people  in  state  (r,t)  being  promoted 
at  time  t to  state  (r+1, t+1) 

1-f [(r,t); (r,t+l)]  - f [(r,t); (r+1, t+1)]  * fraction  of  people  in  state 
(r,t)  being  separated  from  service  at  time  t 
c[(r,t)]  * cost  of  having  a person  serve  his  year  t in  rank  r 

e[(r,t)]  ■ effectiveness  or  value  of  having  a person  serve  his  year 

t in  rank  r. 

We  call  the  above  description  of  the  model  the  state  model.  Figure  1 shows 
the  possible  states  and  transitions  among  them  for  a state  model  in  which  R = 3 
and  T * 5.  After  making  the  assumption  that  effectiveness  as  well  as  costs 
are  additive,  Grinold  goes  on  to  state  the  problem  of  optimal  manpower  system 
design  as  a linear  programming  problem  of  either(a)  maximizing  total  effective- 
ness subject  to  cost,  time  in  rank,  minimum  force  size,  average  length  of 
service,  and  other  constraints;  or  (b)  minimizing  total  cost  subject  to 
minimum  effectiveness,  and  other  constraints  similar  to  the  above.  He  goes 
on  to  suggest  that  sensitivity  analysis  of  the  linear  program  would  supply 
valuable  information. 

In  [2,3]  Chames,  Cooper,  Lewis,  and  Niehaus  show  how  a linear  program- 
ming manpower  system  model  can  be  approximated  by  means  of  capacitated  network 
model.  Here  we  shall  adapt  their  suggestion  to  Grinold's  model. 

The  first  step  is  to  relabel  the  states  using  a lexicographic  ordering. 
Figure  2 shows  the  new  and  old  listing  of  states  for  the  example  of  Figure  1. 

Note  that  the  retirement  and  separation  states  are  numbered  K+l  and  K+2, 
respectively. 

The  network  with  the  new  state  labels  is  shown  in  Figure  3.  For  the  time 
being  ignore  the  numerical  exanple  which  is  also  imposed  on  that  figure. 


-4- 

We  now  introduce  the  notation  necessary  for  the  transshipment  model. 
First  the  "from"  states 

I = {1,2,. ...K}; 

then  the  "to"  states 

J = {2,...,K,K+l,K+2}  . 

We  call  state  1 the  inlection  state;  states  2 through  K are  called  trans- 
shipment states  because  they  appear  both  as  "from"  and  "to"  states;  finally, 
state  K+l  is  the  retirement  state  and  state  K+2  is  the  separation  state. 
Next  we  define 

* manpower  flow  from  state  iel  to  state  jeJ 

These  variables  are  bounded  by 

0 if  manpower  cannot  flow  from  state  i to  j 
upper  bound  on  the  flow  when  flow  is  possible. 

These  upper  bounds  are  varied  to  achieve  the  desired  fractions  of  people 
being  promoted,  continued  in  rank,  separated,  etc.  of  the  state  model.  We 
also  define  initial  and  final  flows 

Xq  « yearly  injection  of  new  employees  into  state  i 

x^  * yearly  employee  retirements 

xg  * yearly  employee  separations. 

For  feasibility  we  assume 

x0  " Xr  + xs  (1) 
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With  the  above  notation  we  can  now  state  the  feasibility  constraints 
for  the  transshipment  model: 


E x * x for  all  iel  (2) 

jeJ  1J 


E 

iel 


'ij 


xQ  for  all  jej-{K+l,K+2} 


(3) 


E 

iel 


Xi,K+l 


XR 


(4) 


E 

iel 


Xi,K+2 


(5) 


0 £ for  all  iel,  jeJ. 


(6) 


If  i is  a transshipment  state,  that  is,  i e I fl  J * {2,...,K}. 
then  the  variables  x are  called  transshipment  variables.  If  index  i cor- 
responds to  state  (r,t)  the  transshipment  variables  have  the  following  special 
interpretations : 


Xq-x^  ■ the  number  of  employees  in  rank  r at  time  t 
x^  * the  number  of  employees  not  in  rank  r at  time  t. 


Obviously  employees  not  in  rank  r at  time  t are  either  in  other 
ranks,  or  have  retired,  or  have  been  separated. 

We  turn  now  to  the  statement  of  objective  functions  for  the  model.  We 
shall  consider  here  just  two  objectives,  minimize  cost,  or  maximize  effectiveness. 
For  the  cost  objective  we  define 


« 
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t 


'ii 


for  i e I f"l  J 


(7) 


'ij 


!«  tf  state  i is  not  connected  by  an  arc  to  state  j,  or 


yearly  cost  of  having  an  employee  transiting  from 
iel  to  jeJ  where  i ^ j 


(8) 


'i ,K+1 


!®  if  retirement  is  impossible  from  i,  or 

!total  retirement  cost  if  retirement  is 
^ possible  from  i 


(9) 


Ci  K+2  = cost  separating  an  employee  from  state  i (10) 

With  these  definitions  the  minimum  cost  transshipment  objective  function 
is 

Minimize  {c  = Z £ ex}.  (11) 

iel  jeJ  1J  iJ 

The  problem  stated  in  ( 2 ) - (6 ) and  (11)  is  a transportation  problem, 
which  is  a linear  programming  problem  having  special  structure.  The  dual 
problem  is  given  by 

Maximize  Z a . u + E b v (12) 

iel  11  jeJ  J J 

Subject  to 

u^  + Vj  - < c^  for  iel  and  jeJ  (13) 

w^j  > 0 for  iel  and  jeJ  (14) 

where  the  dual  variables  u^  for  iel  are  associated  with  constraints  (2),  the 
variables  v^  for  jeJ  are  associated  with  constraints  (3),  (4),  and  (5),  and 

the  variables  w^  for  iel  and  jeJ  are  associated  with  constraints  (6). 


-7- 


For  a given  basis  B it  is  possible  to  determine  a one-parameter 


family  of  solutions  and  v ^ to  the  equations 


ui  + Vj  “ Cij  f°r 

such  that  d . . * u.  + v.  is  unique, 
ij  i j 

of  the  dual  solution  given  by  a basis 
shown  that 


e B (15) 

The  matrix  D = [d.  .]  the  dual  matrix 
ijJ  

E.  As  noted  in  [9]  it  can  easily  be 


* max(0,ui  + vj  " cij)  for  ieI  and  jeJ 


(16) 


for  any  dual  feasible  basis  B.  Therefore  the  dual  variables  u^  and  v^ 
suffice  to  characterize  the  dual  solution. 

For  the  effectiveness  objective  we  define  an  effectiveness  function 
e^j  on  arcs  (i,j)  as  follows: 


e^  * 0 for  i e H J 


(17) 


r = - ® if  state  i is  not  connected  by  an  arc  to  state  J,  or 

» yearly  effectiveness  of  having  an  employee 
transiting  from  state  iel  to  jeJ 


1 ,K+1  * ° £°C 

iel 

(19) 

L , K+2  ' 0 £or 

iel 

• 

(20) 

The  problem  of  determining  the  effectiveness  of  single  Individuals  has  been 
addressed  by  using  Delphi  techniques  by  J.  R.  Schmidt  and  R.  K.  Hovey  in  [7]. 

They  had  difficulty  extending  these  individual  results  to  effectiveness  of 
groups  of  individuals.  Nevertheless,  we  shall  here  make  the  (admittedly  heroic) 
assumption  that  effectiveness  is  an  additive  function.  Thus  the  total  effective- 
ness objective  is 
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Maximize  fE  = E E e x } . (21) 

iel  jeJ  1J  1J 

It  should  be  remarked  that  In  [5]  Grinold  treats  a case  in  which  the  total 
effectiveness  function  is  non-linear.  As  might  be  expected,  the  computational 
problems  in  that  case  are  difficult. 

The  dual  problem  to  the  transportation  problem  defined  by  ( 2 ) - ( 6 ) and 
(21)  is  just  like  the  one  given  by  (12)- (14)  except  that  the  word  "Maximize" 
in  (12)  must  be  replaced  by  the  word  "Minimize."  The  same  remarks  made 
previously  about  the  dual  problem  apply  here. 


3 . AN  EXAMPLE 

A numerical  example  with  cost  minimization  objective  function  is  shown 

in  network  form  in  Figure  3 and  in  transshipment  form  in  Figure  4.  Actual 

numbers  were  chosen  for  c and  for  each  arc;  also  Xq  = 1000, 

x * 681  and  x = 319  were  chosen.  From  the  actual  flows  in  the  solution 
R b 

to  the  problem  it  is  easy  to  calculate  the  fractions  of  people  retained  in  rank 
f [(r,t); (r,t+l)] , the  fractions  of  people  being  promoted  f [(r,t); (r+l,t+l)] , 
and  from  these  two  the  fractions  of  people  being  separated  at  each  state.  The 
upper  bounds  must  be  selected  by  trial  and  error  to  (a)  achieve  feasibility, 

and  (b)  to  achieve  approximately  the  desired  promotion  fractions  in-rank 
fractions  at  each  state.  Charnes , Cooper,  Lewis,  and  Niehaus  selected  upper 
bounds  to  achieve  similar  fractions  in  their  models  [2,3]. 

It  is  obvious  that  with  the  minimum  cost  objective,  the  optimal  solution 
will  try  to  make  the  following  decisions  on  the  arcs  leading  from  a given 


state  i: 
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(a)  Separate  as  many  persons  as  possible;  i.e.,  put  maximum 
flow  possible  on  the  upward  slanting  arrows  in  Figure  3; 

(b)  Retain  as  many  people  in  rank  rather  than  promote  them; 
i.e.,  put  maximum  flow  possible  (after  satisfying  (a))  on 
the  horizonal  arrows  in  Figure  3.  (This  remark  does  not 
hold  on  the  bottom  row  of  the  figure  where  no  promotion 
is  possible.) 

From  these  remarks  it  follows  that  the  basis  cells  of  the  optimal  solution 
to  the  transshipment  problem  in  Figure  3 will  consist  of  (i)  the  transshipment 
cells  and  (ii)  the  promotion  cells,  plus  other  cells  as  necessary  to  complete 
the  basis.  The  other  cells  include  the  in-rank  cells  at  the  highest  rank, 
since  promotion  is  not  possible  from  these  states,  and  the  retirement  cells 
in  column  13,  since  again  no  other  alternative  is  available.  (There  is  one 
additional  cell  (5,14)  needed  to  complete  the  basis.)  These  conclusions  can 
be  confirmed  by  looking  at  the  circled  basis  cells  in  Figure  4 and  the  upper 
bounded  cells  in  that  figure  indicated  by  squares. 

From  these  remarks  it  might  and  can  be  concluded  that  the  solution  to 
the  transportation  problem,  or  the  network  problem  of  Figure  3 is  not  difficult. 
This  is  true, and  advantage  can  be  taken  of  the  fact  that  good  starting  bases  are 
easy  to  find,  in  the  development  of  special  algorithms  for  solving  the  problem. 

With  the  maximum  effectiveness  objective  function,  the  corresponding 
obvious  optimal  decisions  are: 

(c)  Promote  as  many  people  as  possible;  i.e.,  put  the  maximum 
flow  possible  on  the  downward  slanting  arrows  in  Figure  3; 

(d)  Retain  as  many  people  in  rank  as  possible  rather  than  fire 
them;  i.e.,  put  maximum  flow  possible  (after  satisfying  (c)) 
on  the  horizontal  arrows. 

From  these  remarks  most  of  the  optimum  basis  cells  can  be  identified  in  advance 
for  this  case  as  was  done  previously  for  the  minimum  cost  objective  case. 
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4.  DUAL  VARIABLE  ANALYSIS 

The  optimal  dual  matrix  for  the  problem  in  Figure  4 is  shown  in 
Figure  5.  We  use  the  results  of  Srinivasan  and  Thompson  [9]  to  interpret 
the.  relevant  entries  of  the  dual  matrix.  Although  every  entry  in  this  matrix 
has  one  or  more  interpretations,  only  the  entries  in  the  column  13  corresponding 
to  retirement  and  the  entries  in  other  rows  that  have  squares  around  them 
corresponding  to  cells  at  their  upper  bounds,  are  of  interest  here.  Also  it 
turns  out  that  the  optimal  basis  subgraph  of  the  network  of  Figure  1 is  needed 
to  make  such  interpretations.  This  basis  subgraph  is  shown  in  Figure  6. 

Notice  that  none  of  the  horizontal  arcs  (except  those  in  the  bottom  row)  are 
in  the  basis  graph;  all  these  arcs,  which  correspond  to  "retain  in  rank"  are 
used  to  the  fullest  capacity,  i.e.,  are  upper  bounded,  and  hence  cannot  be  in 
the  basis  graph.  Similarly  all  the  upward  slanting  arcs,  which  correspond 
to  separation,  except  for  the  one  from  node  5,  are  upperbounded  and  hence  not 
in  the  basis  graph.  Because  all  the  arcs  corresponding  to  promotions  are 
in  the  basis  this  can  be  called  the  mos t expensive  basis  graph . The  reason 
that  the  optimum  solution  to  the  cost  minimizing  problem  leads  to  the  most 
expensive  basis  graph  is,  of  course,  that  all  the  cheaper  arcs  are  already 
utilized  to  their  fullest,  i.e.,  to  their  upper  bounded  capacity,  in  the 
solution.  Hence  only  the  expensive  arcs  involving  promotions  are  available 
to  be  in  the' basis.  There  is  one  exception  to  this  rule,  namely  the  separation 
arc  from  node  5.  But  that  arc  has  to  be  in  the  basis  to  make  the  basis  graph 
be  connected,  i.e.,  touch  all  the  nodes.  It  is  also  not  upper  bounded  in  order 
that  the  solution  can  achieve  the  desired  number  (319)  of  separations. 

We  can  now  state  a general  rule  for  interpreting  the  numbers  in  the 


dual  matrix  of  Figure  5: 
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Dual  Matrix  Rule.  Construct  the  directed  basis  graph  G of 
an  optimal  solution  to  the  transportation  problem;  find  the 
unique  path  in  G from  node  i to  node  j;  then  the  entry 
d^j  of  the  dual  matrix  is  the  directed  cost  of  the  path  from 
i to  j,  i.e.,  it  is  the  sum  of  the  signed  arc  costs  on  the 
path  with  a plus  sign  on  the  cost  when  the  arrow  on  the  arc 
points  in  the  same  direction  as  the  path,  and  a minus  sign  on 
the  cost  when  the  arrow  points  in  the  opposite  direction  as  the 
path. 

For  instance,  consider  the  15  entry  in  the  upper  left  hand  corner  of  Figure  5; 
it  corresponds  to  i = 1 and  j = 2.  In  Figure  6 the  path  from  node  1 to 
node  2 consists  of  nodes  1,  6,  10,  11,  7,  and  2 in  that  order.  The  sum  of 
the  signed  cost  entries  on  this  path  is 

12  + 16  + 17  - 17  - 13  - 15 

as  indicated  by  the  dual  variable  rule.  This  same  path  is  indicated  by  the 
nodes  connected  by  arcs  in  Figure  4 (the  transshipment  nodes  (i,i)  should  be 
ignored,  since  the  cost  on  such  nodes  is  0.) 

As  a second  example,  consider  the  entry  223  in  row  1 and  column  13  of 
Figure  5.  Applying  the  dual  variable  rule  we  see  that  it  should  be  the  sum 
of  the  costs  on  the  path  from  node  1 to  node  13  in  Figure  -6.  Adding  this  gives 

12  + 16  + 17  + 18  + 160  * 223 


which  verifies  the  rule.  In  the  same  way  it  is  easy  to  check  the  following: 
The  entries  in  column  13  give  the  most  expensive  retirement  paths  from  the 
states  corresponding  to  rows.  In  other  words,  these  numbers  give  upper  bounds 
on  the  total  cost  of  having  an  employee  in  each  of  the  employment  states. 


i 
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In  order  to  interpret  the  entries  of  the  dual  matrix  marked  by  squares 
in  Figure  5,  which  correspond  to  upper  bounded  cells,  we  make  use  of  the  results 
of  Srinivasan  and  Thompson  [9]  on  bound  operators  for  parametric  programming 
for  the  transportation  problem.  Specifically  we  make  use  of  the  definition 
of  a bound  operator  on  p.  221  and  Theorems  9 and  10  on  p.  223  of  that  reference. 
The  results  given  there  can  be  briefly  summarized  as  follows:  Suppose  at 
cell  (p , q ) the  upper  bound  U is  changed  to 

pq 


u' 

pq 


U +6 

pq 


(22) 


where  5 is  a positive  or  negative  number  that  is  "sufficiently  small"  in  a 
sense  to  be  defined  shortly,  see  (25)  below.  Let  B be  the  basis  of  an  optimal 
solution  to  the  original  problem.  Then  there  exist  positive  constants 
and  n,  such  that,  if  6 satisfies 

pq 


pq 


- Li  < 6 < LA 

pq  ~ ~ pq 

then  B is  still  the  optimal  basis  to  the  new  problem  which  is  the  same 


(23) 


except  Upq  is  replaced  by  U^.  Moreover,  if  Z and  z'  are  the  original 
and  new  objective  function  values  then 


Z'  - Z if  (p,q)  i UB 


Z - Z - 6d  if  (p,q)  e UB 

pq 


(24) 


(25) 


where  UB  is  the  set  of  cells  for  which  x ■ U , 

pq  pq 

upper  bounds.  The  determination  of  the  numbers  and  a 

pq  pq 


i.e.,  cells  at  their 

is  not  difficult 

pq  pq 

and  is  explained  in  reference  [9]. 

In  essence  what  (24)  says  is  that  if  (p,q)  is  not  in  the  set  of  upper 
bounded  cells  then  changes  can  be  made  to  U within  limits  given  by  (23) 
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without  changing  the  optimal  solution  or  objective  value  of  the  problem. 

What  (25)  says  is  that  if  (p,q)  is  in  the  set  of  upper  bounded  '■ells,  which 
is  true  for  the  cells  marked  with  squares  in  Figure  5,  then  making  changes 
of  type  (22)  cause  changes  of  type  (25)  provided  6 satisfies  (23).  Since 
the  upper  bounds  U can  be  changed  by  the  decision  maker  in  order  to  achieve 
certain  goals  involving  promotion  fractions,  in-grade  constraints,  separation 
constraints,  etc.,  we  see  that  the  latter  dual  variable  interpretations  are 
very  important.  We  summarize  these  in  the  following  rule:  For  the  upper 
bounded  cells  (p,q)  (marked  with  squares  in  Fig.  5).  changes  of  form 
u'  = U +6  imply  cost  changes  of  form  Z ' * Z - 6d  provided  t 

— p q p q--  *■  * ■ ..  . . a ■■  - ...  . — — p q *-  — ■ — 

sati^ies  the  constraints  - u < 6 < u+  . 

pq  - - 'pq 

As  an  example  consider  the  cell  (1,2)  in  Figure  5;  here  d^  = 15 
and  “ 850  as  shown  in  Figure  3.  If  we  now  replace  IJ^  by 

= 850  +5  we  see  that  the  objective  function  cost  change  is  given  by 
Z'  * z'  - 156.  In  other  words,  increasing  from  850  to  851  by  making 

6 ■ 1 decreases  the  objective  function  cost  by  15;  whereas  decreasing 
from  850  to  849  by  making  6 ■ -1,  increases  the  objective  function  cost 
by  15.  The  way  the  optimal  solution  is  changed  in  each  of  these  cases  can  be 
found  by  tracing  out  the  cycle  shown  in  Figure  4.  It  can  also  be  found  from 

the  graph  of  Figure  6 by  tracing  out  the  path  from  node  1 to  node  2.  We  leave 

the  details  of  the  determinations  of  these  paths  to  the  reader  because  of  lack 
of  space. 

The  dual  variable  analysis  for  the  maximum  ef fectlvr  less  objective 
function  (21)  case  proceeds  in  an  exactly  similar  manner,  but,  of  course, 
with  suitable  changes  in  signs.  Here  the  optimal  basis  graph  will 
correspond  to  a "least  effective"  graph  such  as  the  one  shown  in  Figure  7. 
Details  are  not  given  here  for  lack  of  space. 
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5.  EFFECTIVENESS  VERSUS  COST  TRADEOFF  ANALYSIS 

For  most  purposes  neither  the  minimum  cost  nor  the  maximum  effective- 
ness objective  function  previously  considered  are  necessarily  the  most 
appropriate.  The  usual  situation  is  that  a given  amount  of  money  is  budgeted 
for  personnel  purposes,  and  the  objective  is  to  obtain  the  most  effective 
force  that  costs  the  budgeted  amount.  However,  once  that  problem  is  solved 
additional  questions  are  frequently  asked  such  as:  Given  that  a certain 

level  of  effectiveness  has  been  achieved  with  the  allocated  budget  how  much 

\ 

more  (or  less)  effectiveness  can  be  obtained  with  x percent  more  (or  less) 
money?  Or,  what  will  it  cost  to  achieve  a given  level  of  effectiveness? 

The  best  way  to  answer  these,  and  other  similar  questions  is  to  compute 
the  effectiveness-cost  tradeoff  curve.  For  the  model  being  considered  this 
turns  out  to  be  relatively  easy  if  we  use  the  method  given  by  Srinivasan  and 
Thompson  in  [ ll]  for  solving  multiple  objective  transportation  problems. 
Consider  the  problem  P(6)  with  constraints  given  by  (2)-(6)  and  the  follow- 
ing objective  function 

Minimize  [j  = C-6E]  (26) 

where  C is  the  cost  function  of  (11),  E is  the  effectiveness  function  of 
(21),  and  6 is  a nonnegative  parameter.  When  6=0,  the  problem  is 
the  minimum  cost  problem  previously  considered.  As  6 increases  the  optimum 
solution  to  the  problem  gives  a more  costly,  but  also  more  effective  solution. 
In  Figure  8 the  effectiveness-cost  tradeoff  curve  is  shown.  As  Srinivasan  and 
Thompson  have  shown  in  [ ],  this  curve  is  pareto-optimal , is  non-decreasing, 

and  is  piecewise  linear,  see  also  Geoffrion  [ ].  In  that  reference  they 

exhibit  a very  efficient  parametric  programming  method  for  computing  the 


i 
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trade-off  curve.  The  reader  is  referred  to  that  paper  for  details.  As  noted 
in  Figure  8,  points  above  the  trade-off  curve  are  infeasible,  and  points  below 
are  non-optimal.  A curve  having  these  properties  is  said  to  be  pareto-optimal. 

For  any  given  point  on  the  trade-off  curve  the  optimum  basis  graph  will 
be  a directed  tree  but  with  a less  obvious  structure  than  those  exhibited  in 
Figures  6 and  7.  Hence  it  is  not  as  easy  as  before  to  guess  a good  starting 
solution.  However,  the  dual  matrix  rule  of  Section  4 is  still  valid  and  can 
be  used  to  interpret  the  optimum  solution  in  the  same  way. 

6.  CAREER  ANALYSIS 

The  model  shown  in  Figures  1 and  3 can  be  called  a cohort  model  because 
it  traces  throughout  their  careers  the  flow  of  a group  of  persons  all  of  whom 
enter  the  system  at  the  same  time.  We  can  turn  it  into  an  equilibrium 
steady  state  model  by  considering  five  such  networks  each  starting  at  the 
beginning  of  five  consecutive  years  and  look  across  all  five  for  a single 
year.  More  simply  we  can  just  sum  across  all  the  years  in  Figure  3 to  find  the 
numbers  of  persons  in  each  rank.  These  sums  show  that  of  the  5,000  people 
injected  into  the  employment  system  over  a period  of  5 years  there  are  3046  in 
active  service  in  equilibrium.  They  are  distributed  in  the  various  ranks  as 
indicated  in  Figure  9.  Note  that  there  are  93%  in  rank  1 positions,  6 per'cent 
in  rank  2 positions,  and  1 percent  in  rank  3 positions.  • 

In  the  same  way,  if  we  look  at  the  solution  in  Figure  3 to  find  the 
eventual  fate  of  newly  entering  employees  we  see  from  Figure  10  that  32%  are 
separated  before  retirement,  48  percent  retire  from  rank  1 positions,  17 
percent  from  rank  2 positions,  and  3 percent  from  rank  3 positions.  Also 
shown  in  Figure  10  are  the  conditional  percentages  of  employees  retiring  from 
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each  rank,  given  that  they  are  not  separated  before  retirement.  Note  that 
these  percentages  are  71  percent  from  rank  1,  24  percent  from  rank  2,  and 
5 percent  from  rank  3.  When  compared  to  the  corresponding  numbers  in 
Figure  9,  these  percentages  may  seem  somewhat  paradoxical:  for  instance, 
although  only  1 percent  of  active  employees  in  any  year  are  in  rank  3 
positions,  we  see  that  5 percent  of  those  retiring  in  any  year  do  so  from 
rank  3 positions.  Similar  comparisons  can  be  made  from  each  of  the  other 
ranks.  The  paradox  can  easily  be  explained  by  considering  the  fact  that 
retirees  are  chosen  from  the  oldest  segment  of  current  employees  and  they 
are  more  likely  to  be  in  the  higher  ranks. 

Still  another  use  of  the  cohort  flow  model  in  Figures  1 and  3 is  that 
of  successful  career  analysis.  By  a successful  career  we  shall  mean  one  that 
ends  in  retirement  (rather  than  separation).  The  following  very  simple 
algorithm  can  be  used  to  count  the  number  of  successful  careers. 

(a)  Mark  the  retirement  state  with  a 1. 

(b)  Mark  each  state,  all  of  whose  successors  are  already 
marked,  with  the  sum  of  the  marks  of  its  successors. 

(c)  Stop  when  the  injection  state  is  marked.  Its  mark 
gives  the  number  of  successful  careers  that  are  possible. 

The  calculations  of  the  algorithm  when  applied  to  the  network  in  Figure  1 are 

given  in  Figure  11.  They  show  that  11  successful  career  paths  are  possible 

from  the  injection  state  to  the  retirement  state.  These  successful  career 

paths  are  most  easily  summarized  by  indicating  the  sta  tes  at  which  promotions 

take  place,  as  indicated  in  Figure  12.  The  career  marked  1 means  that  the 

person  entered  employment  at  injection  node  1 and  was  never  thereafter 

promoted.  Career  1,9  Indicates  a career  starting  at  1 going  to  4,  then  to  9, 


i 


and  then  to  retirement.  Etc.. 


In  order  to  get  an  idea  of  how  common  each  of  these  successful  careers 
is  we  trace  each  one  in  Figure  3 and  record  the  arc  along  which  the  fewest 
people  can  travel  in  the  career  because  it  has  the  smallest  optimal  flow  of 
any  arc  on  the  path;  such  arc  capabilities  are  called  bottleneck  flows. 

Each  of  these  bottleneck  flows,  and  its  percentage  of  the  sum  of  all  the 
flows  is  marked  in  Figure  12.  Although  these  percentages  are  not  an  accurate 
evaluation  of  the  difficulty  of  following  a given  career,  they  do  give  some 
indication  of  the  relative  difficulty  of  a given  entrant  following  each  of 
the  11  successful  career  paths.  In  a real  employment  situation,  these  per- 
centages could  be  compared  with  actual  frequencies  of  the  successful  career 
paths.  For  instance  I.  P.  Sherlock  reports  in  [8]  the  proportions  of  people 

in  each  rank  for  certain  trades  in  the  Canadian  armed  forces. 

\ 

7.  CONCLUSIONS 

% The  goal  of  this  paper  was  to  demonstrate  the  usefulness  of  a network 
transshipment  model  to  approximate  a linear  programming  model  for  manpower 
planning  and  design.  The  main  advantages  of  the  network  model  are:  fast 
computing  speeds,  easy  computation  of  effectiveness-cost  trade-off  curve, 
the  availability  of  much  managerial  information  contained  in  the  dual  matrix, 
and  its  usefulness  in  analyzing  career  paths.  The  only  appreciable  disadvan- 
tage is  the  difficulty  af  imposing  extra  constraints.  In  any  specific  instanc 
the  latter  difficulty  can  usually  be  overcome  by  employing  parametric  program- 
ming procedures  as  illustrated  above. 
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Number 

State 

1 

(1,1) 

2 

(1,2) 

3 

(1,3) 

4 

(1,4) 

5 

(1,5) 

6 

(2,2) 

7 

(2,3) 

8 

(2,4) 

9 

(2,5) 

10 

(3,3) 

11 

(3,4) 

K * 12 

(3,5) 

K + 1 = 13 

Retirement 

K + 2 - 14 

Separation 

Figure  2.  Lexicographic  numbering  for  the 
example  of  Figure  1. 


Note  that  x„  = 1000;  x„  = 681,  and  x„  = 319 


minimum  cost  objective  function 


Figure  11.  Application  of  the  algorithm  for  counting 
the  number  of  successful  careers. 

Note  that  11  such  careers  are  possible. 


-28- 


Career  Paths 

Bottleneck 

Flows 

Percentages 

1 

482 

63 

1,9 

67 

9 

1,8 

72 

10 

1.7 

42 

6 

1,6 

25 

3 

1,8,12 

11 

1.5 

1,7,12 

11 

1.5 

1,6,12 

11 

1.5 

1,7,11 

19 

2 

1,6,11 

19 

2 

1,6,10 

4 

.5 

Figure  12.  Career  paths,  bottleneck  flows 
and  bottleneck  flow  percentages 


